Quantum-field dynamics of expanding and contracting Bose-Einstein condensates 
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We analyze the dynamics of quantum statistics in a harmonically trapped Bose-Einstein conden- 
sate, whose two-body interaction strength is controlled via a Feshbach resonance. From an initially 
non-interacting coherent state, the quantum field undergoes Kerr squeezing, which can be quali- 
tatively described with a single mode model. To render the effect experimentally accessible, we 
propose a homodyne scheme, based on two hyperfine components, which converts the quadrature 
squeezing into number squeezing. The scheme is numerically demonstrated using a two-component 
Hartree-Fock-Bogoliubov formalism. 

PACS numbers: 03.75.Nt, 03.75.Mn. 



I. INTRODUCTION 

Multimode quantum fields are the appropriate de- 
scription for a vast array of phenomena in high-energy 
physics, condensed matter physics and cosmology, but 
they are notoriously difficult to analyze theoretically. Un- 
like many quantum fields, it appears feasible that those 
describing degenerate Bose and Fermi gases can be ex- 
perimentally manipulated, detected and studied. The 
development of techniques for this could have broad in- 
fluence on a variety of outstanding problems in physics. 

Understanding the many-body quantum state of a 
Bose-Einstein condensate (BEC) is relevant for sys- 
tems including: squeezed atom- lasers [H, simulations of 
cosmolosdcal particle production in the early universe 
@, H, 0, ||| , and the quantitative description of collapsing 
condensates due to attractive interactions @, 0, [1] ■ 

The quantum field equations describing expanding and 
collapsing condensates are analogues of those describing 
a quantum field on a curved space-time Q. Learning to 
experimentally manipulate and study the analogue sys- 
tem provided by expanding and collapsing condensates 
promises us a new window on processes such as signa- 
ture change that may be relevant for the evolution of the 
early universe Q . The quantum field of expanding and 
contracting condensates is the subject of this paper. 

While quantum-field models exist to approximately de- 
scribe the time-evolution of Gaussian quantum states in 
such non-equilibrium situations [lClllllll2l.ll3lll4llla.ll6l. 
E3) EH El j they require us to specify the initial quantum 
state of the condensate, which is complicated if interac- 
tions are present 0, HH, [22j . 

Due to the controllability of atomic interactions by Fes- 
hbach resonances, quantum field dynamics in a BEC can 
however be examined starting from a non-interacting ini- 
tial state. Then, we assume that the many-body quan- 
tum state is represented by a coherent state; perhaps 
originating from a mixture of all different phases, as in 
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the optical laser [23j. This initial situation was realized 
in experiments on collapsing Bose-Einstein condensates 
with attractive interactions [24j. We also theoretically 
consider it here, but with interactions suddenly rendered 
attractive or repulsive. 

We show that the condensate's state evolves from co- 
herent to quadrature squeezed due to the Kerr effect. 
For repulsive interactions a single-mode model provides 
a qualitative description of the squeezing. However, a 
multi-mode analysis is required for a quantitative de- 
scription, especially in the attractive case. 

The experimental detection of quadrature squeezing 
requires a phase-reference, such as in homodyne detection 
[2a] . We propose such a scheme based on a splitting of the 
Bose-Einstein condensate into two separately conserved 
hyperfine components. We demonstrate the proposal us- 
ing two-component Hartree-Fock-Bogoliubov (HFB) the- 
ory with realistic experimental parameters. The homo- 
dyne detection is found reliable even in the presence of 
imperfections like interactions between the local oscilla- 
tor component and the squeezing field. Also the Kerr 
squeezing of the local oscillator itself does not prevent 
us from attaining a measurable reduction of the number 
variance. 

A single mode model of Kerr squeezing has previously 
been found useful despite the presence of multi-mode 
effects [l|. The authors of Ref. [l[ also observe that 
the interference of two quadrature squeezed atom-lasers 
can yield a number squeezed state, in accordance with 
our successful simulation of the homodyne scheme with 
squeezed local oscillator. In contrast to Ref. [l|, which 
is focussed on the creation of a squeezed atom laser, we 
consider a simpler setup and target studies of quantum 
field dynamics. A core ingredient in our work is a Fes- 
hbach resonance, allowing the use of a simple coherent 
initial state. 

This paper is organized as follows. Section [Til provides 
a brief overview of the three quantum theories employed: 
the single mode model, HFB theory and the truncated 
Wigner approximation (TWA). In section [HT] we analyse 
Kerr squeezing in a harmonically trapped condensate and 
in section IIVI propose a matter- wave homodyne scheme 
to detect it. Technical details regarding variances in the 
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HFB formalism as well as the two-component HFB equa- 
tions of motion can be found in the appendix. 



II. METHODS 

An ensemble of Bose condensed atoms in a harmonic 
trap is described by the many-body Hamiltonian 

,t 



where 



H= /^x* T (x) J ff (x)*(x) 



(1) 



(2) 



is the single particle Hamiltonian, and \P(x) denotes the 
field operator in the Heisenberg picture that annihilates 
atoms of mass m at position x. We have assumed a 
contact interaction of strength Uq — 4Trh 2 a s (t)/m with 
time dependent scattering length a s (t) and a spherically 
symmetric harmonic potential U(x) — tocj 2 x 2 /2. Note 
that we use the physical coupling for the parameter Uo 
directly, rather than the bare coupling, which we justify 
in section III Bl The Heisenberg equation for the field 
operator is: 

ih^f± = H *( X )+U &(x)*(x.)*(x). (3) 
at 

In the following subsections we briefly introduce several 
methods to obtain approximate solutions to this multi- 
mode quantum field problem. 



A. Single mode Kerr squeezing 

Among the single particle bases in which we can an- 
alyze the atom-field dynamics of Eq. ([3]), the harmonic 
oscillator basis takes a special role. The condensate is 
initially assumed to be non-interacting and in the trap 
ground state. Thus in the oscillatory basis we assume it 
is in a coherent state of one single particle mode. 

Let us expand the field operator as: ^(x) = 
J2kLo ( ^fe( x )ofe- The ^fc(x) are the eigenstates of the sin- 
gle particle Hamiltonian for a harmonic potential, with 
a collective index k labelling all quantum numbers. The 
operator a,k annihilates an atom in eigenstate (^(x), with 
-ffc)C/9fc(x) — hu>k l fik('X-)- Using this expansion, we can 
rewrite Eq. ([3]) as 



th-a k 
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The Ukimn are overlap integrals of the form 



Uklr 



U /d 3 x^(x)^(x)^ m (x)^„(x) 



(4) 



(5) 



Initially all the atoms are in state ^o( x ), and for short 
times we can approximate Eq. ([¥]) by: 



S 9 - 



hwo&o + Uooooala aQ 



(6) 



Using the trap ground state ^o(x) = Aexp [— r 2 /(2a 2 )] 
with A = (ttct 2 )- 3 / 4 , er = (hm/uu)- 1 / 2 , we find U 00 oo = 
[/o(27rer 2 )~ 3 / 2 . The ground state energy term in Eq. ([6]) 
can be eliminated using rotating frame operators: b = 
exp (ibjQt)a,Q. The equation of motion then becomes: 
i-§lb = — x&bb, with \ = Uoooo/h. It is known that 
this operator equation gives rise to Kerr squeezing [25j 
in the evolution of the quantum state. Further details 
can be found in Ref. Here we merely state the most 
important facts. 

One can define quadratures for the state </?o( x ) as: 



X 



a e 



-if) 



(7) 



The variance of these operators, [AX 8 ] 2 = ([X 9 ] 2 ) - 
(X s ) 2 , gives information about the shape of the quan- 
tum state's Wigner functon in phase-space [25[. In the 
Kerr effect the variance in one quadrature 9 sqz drops be- 
low the value for a coherent state ([AA e ] 2 = 1), while 
for the orthogonal quadrature 9 sqz + ir/2 it increases. 



B. Hartree-Fock-Bogoliubov theory 

To go beyond the previous section and include multi- 
mode aspects of the quantum evolution, we make use 
of the HFB formalism [TFj, HH, E3, El- Thus we de- 
compose ^(x) into a condensate part </>(x) and quan- 
tum fluctuations x(x), such that vE'(x) = </>(x) + x(x) 
and (^(x)) = 0(x). The quantum fluctuations can be 
described in terms of their lowest order correlation func- 
tions: the normal density Gat(x, x') = (x^(x')x(x)) and 
anomalous density Ga(x,x') — (x(x')x(x)). The result- 
ing equations of motion and their implementation for a 
spherically symmetric, trapped condensate have been de- 
scribed in Refs. 

We now explain how to calculate variances of the oscil- 
lator ground state quadratures ([7]) in HFB theory. Since 
we have dk = Jd 3 x ( y 9^(x)*I'(x), we can use 

(a\ay) = /d 3 x fd 3 x>^(x')^(x)(ft(x')$(x)) (8) 



Expressions like Eq. ([8]) are all we need to extract the 
quadrature variance of the trap ground state [AA e ] 2 
from our HFB simulations using spherical symmetry (de- 
scribed in [t}). The result for the variance of the 9 
quadrature is: 



[AX» 



6]2 



1 + 2 /d 3 x /d 3 x' 



<A)(x'Vo(x)*Gat(x,x') 



2SRc{^(x')^(x)G,(x,x')e- M } 



(9) 
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For our analysis of the homodyne scheme we are also 
interested in the variance of the total atom number. Us- 
ing N = Je? 3 x ^(x)^(x) we derive: 



d J x /d 3 x' 



[AN]' = Ntot 
2fHe{r(x)0*(x')G A (x,x')} + 20*(x)0(x')G JV (x,x') 



|G A (x,x')| 2 + |G j v(x,x')r 



(10) 



Further details about the form of Eqs. ([9]) and fTO]) that 
we use for numerical solutions in spherically symmetric 
situations, are given in appendix lAl 

In this work we have found that the results are inde- 
pendent of the numerical cutoff K = Tt/Ax, where Ax 
is the grid spacing, only without coupling renormaliza- 
tion. As previously noted [12, HH, the diagonal part 
of G a is the only divergent quantity in the formalism. 
However, this contributes negligibly to the variances of 
interest here, Eqs. © and (JQJ) . 



C. Truncated Wigner theory 




20 30 
t[ms] 



FIG. 1: (Color online) Evolution of the condensate density 
and ground state population for scenarios I (a,c) and II (b,d), 
as described in the text. The ground state population in 
panels (c) and (d) is defined by No = Jd 3 x ipo(x)*<^)(x). In 
either case the mean field undergoes visible change of shape 
but nonetheless most of the population remains in the ground 
state (c,d). Note the different time-scales for the two scenar- 
ios, and that the spatial axis is a radial coordinate. 



We have highlighted the value of verifying numerical 
quantum-field theory results for a BEC by using two 
quite different formalisms in Ref. 8] . Here we follow the 
same approach, by investigating quadrature squeezing 
using the HFB method as well as the truncated Wigner 
approximation (TWA) 0, M, ES EE d. We havc 
given a compact summary of the method and its imple- 
mentation in Ref. [8J. Our TWA simulations are based 
on solutions of the relevant stochastic differential equa- 
tion in the harmonic oscillator basis 18]. The theory is 
then expressed in terms of the stochastic wave function 

a ( x ) = E„ a ™VnW. 

To determine the quadrature variances in this frame- 
work, we use the appropriate correspondences between 
averages of the stochastic wavefuntions and operator ex- 
pectation values. Most importantly [ToT ]: 

<^ -> \ (i^Jfm) + (*mH)) (11) 

Using Eq. (fTTj) we obtain 

[Al 9 ] 2 =2 (^- |arJ| 2 ) 

+ 2me{(a^-o^ 2 )e- 2ie } . (12) 

III. SQUEEZING OF THE GROUND-STATE 

In this section we discuss our simulations of ground 
state quadrature squeezing in a harmonically trapped 
condensate using the TWA and HFB. We begin with a 
87 Rb condensate of 6000 atoms in a spherical trap with 
u> = 12.8 x 27T Hz. Initially, the interaction between the 



atoms is tuned to zero by use of a Feshbach resonance 
so that all atoms occupy the trap ground state. The 
scattering length is then suddenly switched to an cither 
attractive or repulsive value, which we label ad yn - 

We consider two cases: scenario I with Qdyn — -12a 
and scenario II with ad yn = +12ao. In scenario I the con- 
densate contracts, as expected. We evolve it for 5 ms, 
which is less than the Gross-Pitaevskii/HFB collapse 
time [24| of about 7.5 ms for this interaction strength 
0. More than 96% of the population stays in the trap 
ground state mode for these 5 ms. All this is shown in 
Fig. [TJ For the repulsive interactions in scenario II, the 
BEC is stable and we evolve it for the longer span of 
45 ms. During this time the cloud performs roughly one 
breathing oscillation as shown in Fig. [1] (b). Less popu- 
lation is transferred to non-ground state modes than in 
the attractive case. 

We determine the evolution of the quantum state of 
the condensate in the HFB approximation. The initial 
state is a coherent state, with Ga = Gn = 0, and 
these correlation functions evolve nonzero values describ- 
ing the Kerr squeezing. Our numerical results are shown 
in Fig. [21 They show significant quadrature squeezing, 
which evolves to a maximum and then decreases. We 
find the largest squeezing in the repulsive case (Scenario 
II): up to 9 dB. These results are compared with the pre- 
dictions of the single mode model of section III Al For 
the repulsive scenario II the single mode model approxi- 
mates the maximal squeezing attained and the time scale 
on which it occurs. For the attractive scenario I the sin- 
gle mode model is accurate only for short times. This 
is due to the contractive dynamical instability resulting 



in greater production of uncondensed atoms than in the 
repulsive case [1, H3] . 



IV. MATTER- WAVE HOMODYNE SCHEME 

Experimentally one can measure the uncertainty of 
observables, such as the total atom number in a BEC, 
by determining the shot-to-shot variance. However, the 
Hamiltonian ([1]) commutes with the number operator 
N = fd 3 x so it, and its statistics, are conserved. 

The squeezing described in section Hill must therefore oc- 
cur along some quadrature angle 6 ^ of Eq. (0 and 
cannot be experimentally measured without the phase- 
reference provided by the local oscillator of a homodyne 
scheme j25[. Homodyne detection is well established in 
quantum optics. Its core ingredient is a strong laser beam 
whose quantum state is coherent. Its phase provides 
the reference necessary to extract the noise amplitude 
in any quadrature. In this section we propose a homo- 
dyne scheme for Bose-Einstein condensates in a harmonic 
trap, making use of interference between different hyper- 
fine states of the condensed atoms. 

To this end, we consider a condensate with two com- 
ponents denoted by |1) and |2). Atoms are converted 
between components by applying electromagnetic fields; 
for example using microwave and RF fields in the 87 Rb 
experiments of Ref. [28|, [29[. The Hamiltonian for our 
two component Bose gas is hence 

J i=l,2 v ' 

+ Yl + ^I^ + fT^i}. (13) 

The field operator ^(x) carries a hyperfine index i = 1, 2 
and its spatial argument is suppressed in Eq. (|13[) . We 
assume identical traps for both components. There are 
two intraspecies (J7n, U22) and one interspecies {U\i = 
U21) interaction strength. The magnitude of the coupling 
coefficient |fi|, the Rabi frequency, controls the rate of 
conversion between the species. The phase of fi controls 
the relative phase of the two atomic components (see Sec. 

EES). 



A. Homodyne detection 

Let i& denote an atomic field, represented by compo- 
nent |1), that undergoes Kerr squeezing. The spatial ar- 
gument x is suppressed and we will denote \J/(x') by 
Let $ denote the local oscillator atomic field, represented 
by component |2), which we assume to be at all times in 
a coherent state ($(x)) = b{x)e lS ^ x \ with 6(x) and 0(x) 
some real functions. We also assume a large amplitude 




t [ms] t [ms] 



FIG. 2: Minimum quadrature variances V = mine [AX 9 ] 2 
corresponding to maximal squeezing versus evolution time. 
We compare HFB simulations (x, Eq. (lA4jH with truncated 
Wigner results (solid line, Eq. (|12p . the dotted lines indicate 
the sampling error), (a) Scenario I. (b) Scenario II. Both 
panels also include the corresponding analytical result for the 
single mode Kerr effect (dashed). We use Eq. (1) of Ref. 0], 
with x = 0.0134. 



for the local oscillator, 

\b(x)\ 2 » ($t(x)$(x)>. (14) 

In the following we denote X 6 = ^>e~ 10 + SSf'e , where 
0(x) is abbreviated to 6, suppressing the spatial depen- 
dence, and 0(x') is abbreviated to 6'. 

We denote the combined field (p = *P + and deter- 
mine its number variance; [AAtot] 2 = (N 2 ot ) — (N tot ) 2 for 
Ntot = jd 3 ~x. <fi ( x )v?( x )- The total number uncertainty 
can be obtained from: 

(N 2 ot ) = jd 3 x jd 3 x' 

x ((*+$y (v'+&y 

= Jd 3 x Jd 3 x' {(^Mt'*') + (¥)b 2 'be l9 
+ {¥')b 2 b'e l8 ' + (i>)b 2 'be- ie + {^')b 2 b' e- ie ' 

+ b 2 {¥'i>') +6 2 '(* t #) +bb'e- l(9 '- 9) {i>^ 1 ) 

+ bb'e-W- 8 ') (M?) + bb'e-^+^^if') 
+ bb'e i{ - 9+e '\¥'¥) + b 2 b 2 ' + 6e- ifl 

+ b'e ie '{&m r ), (15) 

and 

(JV tot ) = Jd 3 x (<*t$) + b( x )(X e ) + 6(x) 2 ) . (16) 

In writing Eq. (fT5|) we have factored the correlations be- 
tween the fields ^ and <K With a strong local oscillator, 
Eq. (|14p , we need only retain the leading order in 6, and 
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obtain: 



These have the solutions I3C 



[AiV to t] 2 = Jd 3 x yd 3 x'6(x)6(x') 

x { e -<(»-«') [<5 3 (x-x') + Gat(x,x')] + e^-^G^x', x) 
+ e- W)G A (x, x') + e^+^G^x, x')} 



*i = *i(0) cos(|n|t/fi) - «p * 2 (0) 8in(|fi|t/ft), 



= iV to t(l + 2 yd 3 x /d 3 x' ^ (x)^o(x') 
■Xtljtfx ^x' e- i; ( 9 ' +e )G A (x,x' 



-t(e'-e) 



Gjv(x,x') 



(17) 



For the last equality, we have assumed that the local 
oscillator is in the trap ground state, so that &(x) = 
ipo(x)\/ N 2 — 9?o( x )v / -^tot- N2 is the atom number in the 
local oscillator and we have noted that the local oscilla- 
tor is highly populated compared to the squeezed field. 
If we also assume that the phase of the local oscillator is 
homogenous 0(x) = 0(x') = 8 we can write: 



[AN tot ] 2 = N tot [AX^ 



(18) 



using Eq. ©. The quadrature angle 9 is here given by the 
phase angle of the local oscillator condensate. Through 
the mixing of the local oscillator with the squeezed field, 
the quadrature of reduced uncertainty can be rotated into 
the number "quadrature" . 



B. Component mixing 

The experimental scheme we are modeling uses two 
pulses of radiation resulting in a non-zero fl in Eq. (|13|) . 
The first creates a small component |1) field from an ini- 
tial, large, non-interacting component |2) field. Compo- 
nent |1) is the field whose squeezing we seek to measure, 
and component |2) serves as the local oscillator. Between 
the first and second pulses the Kerr squeezing evolves, 
due to self-interaction in the |1) component. The second 
pulse, a time t ovo i vc after the first, mixes the target field 
and the local oscillator so that the quadrature squeezing 
may be inferred from the number variance. 

The particular quadrature measured is determined by 
the phase of the complex coupling O. Since this is the 
phase of an RF field, it may be easily adjusted between 
the two pulses, allowing access to all quadratures. 

We now consider the electromagnetic coupling in more 
detail. When it is switched on, the atoms undergo Rabi 
oscillations between the hyperfine components. These 
can be understood by isolating the relevant parts of the 
Heisenberg equations for the field operators 



* 2 = *a(Q) oos(|fi|t/R) 



101 ~ 

i^*i(0)sin(|O|t/ft). (20) 

After t = = 7wr/(4f2) the trigonometric functions 
have the value 1/V2 and each field has equal contribu- 
tions from the initial fields, called a ir/2 pulse, with a 
relative phase determined by the phase of the coupling 
coefficient fl. The modulus of fi is chosen such that 



C 7T/2 



< t, 



evolve • 



Feshbach resonances 



(19) 



The situation described in section llV Al will be difficult 
to achieve. Since we require more atoms in the local os- 
cillator component |2) than in the component in which 
we wish to measure Kerr squeezing |1), this configura- 
tion will give rise to Kerr squeezing in the local oscillator 
itself due to its self-interaction (t/22), invalidating the as- 
sumption for it to be in a coherent state. Further, there 
are interactions between the local oscillator and the com- 
ponent to be measured (U12), which can even result in 
spatial phase separation [31| , and will also affect compo- 
nent one's Kerr squeezing. 

Ideally, we would like to set U22 — U12 — using 
a Feshbach resonance. However, this would require a 
rare coincidence of resonances in two different scattering 
channels. Hence we also consider three other options for 
improving the situation with a Feshbach resonance: (i) 
increasing Un, possibly to the point that U22 and U12 
can be neglected in comparison, (ii) turning off U12 only, 
and (iii) turning off U22 only. 

We will present detailed results for option (i) and 
briefly comment on our findings for the other options, 
including the ideal case. 



D. Two-component Hartee-Fock-Bogoliubov 
theory 

To investigate the ideas of the previous sections in a full 
multi-mode quantum field simulation of the homodyne 
scheme, we employ HFB theory. For the two-component 
case the simulation uses two condensates, <f>i — (*»), and 
six correlation functions, as detailed in appendix[B] That 
appendix also gives our HFB equations of motion. 

We simulate the evolution of squeezing in component 
|1) followed by mixing with a highly populated compo- 
nent 1 2). We show that this reduces the total number 
variance in one component below the shot-noise limit, 
given by the number in that component. 

Since bothA>i = /d 3 x#|*i andiV 2 = Jd 3 x*^ 2 

com- 
mute with the Hamiltonian (fT3|) for O = 0, we expect the 
number variance of each component, Eq. (|10j) . to change 
only during the mixing step. 
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FIG. 3: (a) Initial shape of "squeezing" condensate |0i(r)| 2 
(solid) and local oscillator 1 02(f) 1 2 (dashed) and the corre- 
sponding (equal) densities after the mixing at t — 0.4 ms 
(dash-dotted), (b) Number variance in condensate compo- 
nent 1 after the mixing step. For arg(f2) = 0.66-7T the Kerr 
squeezing in component one with ming[Xo] 2 = 0.86 is largely 
converted into detectable number squeezing. The thin line in- 
dicates the value of 0.86 for full conversion of the quadrature 
squeezing. 



E. Numerical results 

Now we present two-component HFB simulations of 
the homodyne scheme proposed in section HV A[ We as- 
sume the BEC is initially split into a small cloud in com- 
ponent |1), which is to be squeezed, and a larger cloud in 
component |2), to serve as the local oscillator. We envis- 
age the following creation sequence. After condensation 
the BEC is adiabatically brought to a non-interacting 
initial state with all atoms in the same hyperfine com- 
ponent as done in Ref. ,24]. Using electromagnetically 
induced component mixing, the condensate is then split 
into small and large condensates in different states. This 
is the starting point of our simulations. 

As discussed in section HV C[ after the splitting we as- 
sume that only one of the three couplings [Zy is tuned 
using a Feshbach resonance. For the other two, we used 
the scattering lengths of the 87 Rb hyperfine components 
II) =JF = l,m F = -1) and |2) = \F = 2,m F = I) 
|3ll |32| . which are an = 100. 4ao, 022 = 95.47ao and 
(22i = 98.10cto, where ao is the Bohr radius. We also 
considered parameters appropriate for 85 Rb. For this 
case multi-component scattering length data is not avail- 
able, thus we simply set both non-manipulated scatter- 
ing lengths to a = — 443oo [331 ] . The initial numbers are 
Ni = 600, N 2 = 5400 both in a trap ground state. These 
states are shown in Fig. [3] (a). 

For the results shown in Fig. [5J we increased an by a 
factor of five from its natural value, as in option (i). This 
is not enough to really make U22 and U12 negligible, how- 
ever we found stronger interactions were not numerically 
tractable. In such cases we find that the local oscilla- 
tor becomes squeezed as well, due to its self interaction, 
however \U\\\ > IC/22I ensures that its squeezing does not 
much exceed that of component one. 

For a time tevolve = 0.4 ms we see that the ground- 
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FIG. 4: Quantum state evolution of a two component BEC. 
Panels (a,c) show component |1), (b,d) component 2). The 
variance of the trap ground-state for the quadrature with 
maximal squeezing, Eq. ([9jl, is shown in panels (a,b). The 
variance of the total number in each component, Eq. IjlOj l. is 
shown in panels (c,d) as a fraction of the total atom num- 
ber in each component. For the first 0.4 ms, Q = 0, fol- 
lowed by a pulse of length / x i pu i so = / x 5 fis with 
fl = ?t7r/(4tp U i se ), where / < 1 is adjusted to achieve full mix- 
ing. Subsequently, we again have fl = 0. The two curves use 
a different reference phase: (solid) arg(f2) = 0.66-7T, (dashed) 
arg(Sl) = 0.167T. The grayed region shows the value achieved 
after the component mixing. We did not evolve the system 
past ievolvc + ^pulsc = 0.405 ms, but continued the graph as- 
suming constant variances until 0.6 ms as a visualization aid. 



state of component one, with ad yn = 502ao, develops 
a minimum quadrature variance of about 0.86 (0.66 dB 
squeezing). During this time, since the electromagnetic 
coupling is off, the relative number variances [ANi] 2 /Ni 
remain one to within the Gaussian approximation of the 
HFB method. The slight reduction in the relative num- 
ber variances seen in Fig. |4jc,d) is presumed to be due 
to the development of higher order correlations than can 
not be treated with the HFB, or TWA, method. After 
Evolve, we apply a nonzero coupling for t pu i so ~ 5 /is, un- 
til the components have mixed to equal populations, and 
the total number variance now reflects the quadrature 
variance of the ground state of component one before the 
pulse, for a quadrature angle 9 which is controlled by 
the phase of the coupling, arg(fi). Changing this phase 
we can pick up anti-squeezed or squeezed quadratures, 
shown in Fig. g] (c) and (d) and Fig. [3] (b). 

We find that the homodyne scheme reduces the num- 
ber fluctuations of the atom field after recombination be- 
low N despite squeezing of the local oscillator. Similar 
findings have been reported in [l|. 

In the case shown in Fig. [H corresponding to option 
(i) section HVCi we find that the squeezing does not in- 
crease much with further evolution. Although, it should 
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increase with higher values of an the results presented 
here are sufficient to prove the principle of the homo- 
dyne scheme. For the other two options, (ii) and (iii) of 
section HV C[ we found that setting U22 to zero without 
increasing U\\ results in an even earlier turn-around of 
the squeezing in component one, which therefore is neg- 
ligible. Setting U12 to zero prevents this, but without 
increasing U\\ the local oscillator is much more squeezed 
than the other component. 

Finally, we examined a 85 Rb type scenario, with all 
three scattering lengths negative. For our simulations 
we increased the magnitude of an by a further factor 
of 5. This scenario is feasible if the evolution time is 
much shorter than the collapse time, but shows an earlier 
turnaround of squeezing than the corresponding repulsive 
case. 



V. CONCLUSIONS 

We have shown that harmonically trapped condensates 
with a Feshbach resonance provide a clean and stable 
system to study the dynamics of quadrature squeezing 
in atom-optics. We showed how to implement a matter- 
wave homodyne scheme and numerically demonstrated it 
for experimentally feasible parameters. 

The insight into the squeezing evolution of vacuum- 
fluctuations that our scheme affords might be useful for 
analogue cosmology along the lines proposed in [H, 0] • 

The Kerr effect studied here eventually gives rise to 
strongly non-Gaussian quantum fluctuations [251 ]. These 
have been conjectured to cause a notable discrepancy 
between experiment and theory in the collapse time of 
attractive BECs [8(. This aspect of Kerr squeezing in 
attractive condensates might merit further study. 
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APPENDIX A: QUANTUM VARIANCES IN THE 
HARTREE-FOCK-BOGOLIUBOV FORMALISM 



Here we provide details regarding Eq. ([9]) for a spher- 
ically symmetric situation. We decompose the oscillator 
eigenstates into a radial and an angular part v'nim(x) = 
fni(r)Yi m (9, </>), where Yi m (9,<f)) is a spherical harmonic. 



We also expanded the collective index k — > nlm. Fur- 
ther we choose the spherical polar co-ordinate system r' , 
</>', 9' for the vector x' such that its z-axis points along 
the vector x. 9' then denotes the angle between x and 
x', which appears as argument in the atom field corre- 
lation functions for the case of spherical symmetry 0]. 
Finally we use an expansion of Gjv in terms of Legendre 
polynomials P s : 

M 

G JV (x,x / ) =5^Gr#(r,r')P.(cos0'), (Al) 
where r = |x|, r' = |x'|. With this Eq. ((5|) becomes: 



oc 



("'nlm^n'l'm') = / drr 2 I dr'r 12 

Jo Jo Jo 

M 



00 



2 71 



dcos9 



x 



/ #' / dcos6'Y^Y Vml {6',<l>')YC m {e,(l))P s {cos9') 
Jo J-i s=0 



xG%\r,r')r nl {r)f n , v {r'). 



(A2) 



For the oscillator ground state Y 00 = 1/ \T&K- This allows 
us to carry out the angular integrations and obtain: 

/>oo />oo 

(aS 00 «ooo) = ^ \ drr 2 dr'r' 2 Gf (r,r')f* Q (r)f Q0 (r'). 
Jo Jo 

(A3) 

We have used: J^_ 1 dxP m (x) = 28 m fl. Defining foo{ r ) = 



foo{r)/r and G^,(r, r') = G ( ™/ N (r, r')/rr' we finally ob- 
tain: 



i(m) 



mine 



AX 



= l + 8ir( dr dr'[G^(r,r')/ * (r)/oo(r)] 
v Jo Jo 



dr / dr'[G^(r,r')foo(r)foo(r)} 
Jo 



(A4) 



Inserting the expansion in Legendre polynomials of the 
correlation functions into Eq. (|10p we can also obtain: 



[AAf = N t , 
+ 8tt 2 I dR / dR 1 



lle{<J>*{R)4>*(R')Gf(R, R')} 



4>*(R)4>(R')G^(R,R'] 



G A (R, R') + Gjy (R, R') 



(A5) 



APPENDIX B: TWO-COMPONENT 
HARTREE-FOCK-BOGOLIUBOV EQUATIONS 



To follow the quantum evolution of our two-component 
system through the initial squeezing stage past the mix- 
ing step, we make use of the HFB formalism. The set of 
variables used in 0] must be extended to accommodate 
two hyperfme states. We split the field operators into 
mean and fluctuations: Vti = 0i + xi> ^2 = 02 + X2, 
with (if? n ) = (f>n- We then consider two condensates <jf>i, 
4>2 and six correlation functions (using i £ {1, 2}): 



G J vi(x 1 x')=(xt(x')x i (x)) ) (Bl) 

GAi(x,x') = (xi(x / )Xi(x)) ) (B2) 

Gciv(x,x / ) = (^(x')Xi(x)), (B3) 

G C a(x,x') = <X2(x')xi(x)). (B4) 



^G A1 (x,x') = (iT i + #oi)Givi(x,x / ) 

+ C/ 11 {2[n 1 +ni]G A i(x,x / ) 

+ 7]iGati(x,x')* + 77iGati(x,x') 

+ ?7 1( 5( 3 )(x-x')} (B7) 

+ U 12 {(n 2 +n' 2 )G A1 (x ) x') 

+ CG CA (x',x) + C'G CA (x,x') 

+ £G CT v(x',x)+£'G CT v(x,x')} 

+fi[Gc^(x,x')+G C A(x',x)], (B8) 
Throughout we have used: 



The equations of motion for the condensates and the 
correlation functions Ga/na/ 2 are partially identical to 
those previously presented pj. However, additional terms 
exist in all of them due to the coupling between the 
two hyperfme components. In the following we use the 
notation of Ref. p}, in particular: Gjvi = Gjvi(aO = 
Gati(x,x), 4> n = 0„(x), 4>' n = 0„(x') etc. We fur- 
ther introduce the abbreviations: rij = \(f>j\ 2 + Gjyj, 
Vj = <t>j + Gaj, £, = 4>2<Pi + G C a, C = 0201 + g cn and 
H 0j = -h 2 Vl/(2m) + V(x) + iy J , H oj = -h 2 Vl,/(2m) + 
V(x') + Vj : where we allowed possibly different detun- 
ings Vj (not used in the present work). Our equations of 
motion are then: 



ih4> x =#oi0i + U u ([2G m + |0i | 2 ] 0i + G4i0t) 

+ U 12 (n 2 <pi + G C jv02 + G C a<P 2 ) +^02, (B5) 



(xt(x')x 2 (x)) = G C Ar(x',x)*, (B9) 
(x t 1 (x')x|(x)) = G CA (x',x)*. (BIO) 



One can deduce the equations for component two from 
those of component one by using the following symmetry 
relations, under exchange of particle labels 1^2: 



<-> Gcat(x,x'), 
<-> Gca(x,x'), 



Gcat(x',x)' 
G Cj 4(x',x 

Gcn <-> G CN , 
Gca <-> Gca, 



(B1I) 
(B12) 
(B13) 
(B14) 
(B15) 



The equations of motion for the cross correlation func- 
tions are: 



iHG m fax!) = {Hoi - H' 01 ) Givi(x,x') 

+ C/ u {2[n 1 -n' 1 ]G Jvl (x,x') 

+ ^iGai(x,x')* - t/^Gai^x') j 

+ Ui 2 [ (n 2 - n' 2 ) Gjvi(x,x') 

+ CG C jv(x , ,x)*-C'*G CJV (x,x / ) 

+ ^G C a(x',x)*-^Gca(x,x')} 

+ r!GcAr(x',x)*-fi*G C A'(x,x'), (B6) 



«/lGc7v(x, x') = 

(floi - #02) Gcat(x,x') + 2 [tfimi - EW 2 ] G C jv(x,x') 
+ [/n^G^C^x') - C/ 22 ?7 2 *G C a(x,x') 

+ ^i 2 {(n 2 -n' 1 )G CA r(x,x') 

+ £G A2 (x,x / )*-£'*G A i(x,x') 

+ CG JV2 (x,x')-C'G A ri(x,x')} 

+ n(G JV 2(x,x / )-G J vi(x,x / )), (B16) 



9 



and 

(Hoi + H' Q2 ) G CA (x, x') + 2 [Uum + U 22 n' 2 ] G C a(x, x') 
+ U 11 i] 1 G* CN (x, x') + U 2 2V2 G cn(^, x') 

+ ?7i2-f (n 2 + «i) Gca(x,x') 

+ £ [g^x')* +<5 (3) (x-x')] +e'G A r 1 (x,x') 

+ CG j42 (x,x')+C'*G A1 (x,x')} 



0G A2 (x,x / )+0*G 4 i(x,x / ). 



(B17) 
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